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Abstract. We study the continuum space-time limit of a periodic one dimensional 
array of deterministic logistic maps coupled diffusively. First, we analyse this system in 
connection with a stochastic one dimensional Kardar-Parisi-Zhang (KPZ) equation for 
confined surface fluctuations. We compare the large-scale and long-time behaviour of 
space-time correlations in both systems. The dynamic structure factor of the coupled 
map lattice (CML) of logistic units in its deep chaotic regime and the usual d = 1 
KPZ equation have a similar temporal stretched exponential relaxation. Conversely, 
the spatial scaling and, in particular, the size dependence are very different due 
to the intrinsic confinement of the fluctuations in the CML. We discuss the range 
of values of the non-linear parameter in the logistic map elements and the elastic 
coefficient coupling neighbours on the ring for which the connection with the KPZ-likc 
equation holds. In the same spirit, we derive a continuum partial differential equation 
governing the evolution of the Lyapunov vector and we confirm that its space-time 
behaviour becomes the one of KPZ. Finally, we briefly discuss the interpretation of the 
continuum limit of the CML as a Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) 
non-linear diffusion equation with an additional KPZ non-linearity and the possibility 
of developing travelling wave configurations. 



1. Introduction 

Extended dynamical systems can be modelled with infinite dimensional non-linear 
partial differential equations but also with coupled map lattices (CMLs). The relation 
between the two and, in particular, their behaviour in the spatio-temporal chaotic regime 
is a subject of great interest pQ . 

Whether local chaos can be modelled by stochastic noise is a very deep question 
that has been addressed from different angles and with a variety of analytical and 
numerical approaches. Back in the 40's Von Neumann and Ulam proposed to use the 
logistic map in its chaotic domain as a pseudo random number generator Except 
for a set of initial values of measure zero, the iterates of the logistic map are distributed 
according to p(x) = [7r 2 x(l — x)]~ 1,/2 and the sequence T~ 1 (x) with T(x) = sin 2 (7r/2a;) 
is equidistributed within the unit interval [0,1]. However, it was later realized that 
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the time series thus generated has a number of defects, such as time-correlations and a 
trivial first return-map; the logistic map was subsequently neglected as a pseudo-random 
number generator [31 Ej . 

The question above was later revisited in the context of extended dynamical systems 
modelled with either discrete models or continuous partial differential equations. CMLs 
with short-range spatial interactions and chaotic dynamics of the independent units 
develop collective behaviour and long-range order in finite regions of the phase diagram 
characterised by the strength of the non-linearity and the coupling between the dynamic 
units 5j. Some statistical mechanics notions were shown to be useful to describe their 
spatiotemporal behaviour [fi]-[T3]. 

On the continuum side, in the authors compared the deterministic, 

driven by inherent instabilities, Kuramoto-Sivashinsky (KS) partial differential equation 
equation [2*2*1 "2*""] to the stochastic Kardar-Parisi-Zhang (KPZ) partial differential 
equation [2I]-[2E| for surface growth. A careful study of the scaling properties of 
the solution of both equations showed that they are identical in d — 1 but differ in 
d>2 "T""][""7|, a result confirmed numerically in |19j . 

In this paper we investigate the spatio-temporal behaviour of a one dimensional 
ring of coupled logistic maps by taking a formal continuum space-time limit. First, 
we identify two non-linear terms as leading to effective noise and confinement and we 
argue that the CML has a long wave-length, long time behaviour that is very similar 
to the one of a modified KPZ equation describing surface fluctuations in a confining 
potential. Second, we revisit the equation for the evolution of the Lyapunov vector 
using similar arguments [2Z1- Finally, we reckon that the resulting partial differential 
equation is an extension of the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) [28] 
non-linear diffusion equation with an additional KPZ-like non-linearity and we exhibit 
some travelling wave solutions. 

The article is organized as follows. In Section [21 we recall the definition of the 
logistic map and the CML model we study. In Section El we review the partial differential 
equations we are concerned with. We first describe the KS and KPZ partial differential 
equations and we recall the connection between them in d = 1. We then present the 
FKPP equation and some properties of its travelling wave solutions. In these first 
sections we also set the notation and we explain our goal. In Section [""] we define the 
continuum limit and we discuss each term in the resulting partial differential equation. 
In Section El we present the results of the numerical integration of the CML. We test the 
'noise' term, we confront the behaviour of several observables to the corresponding ones 
in KPZ one dimensional growth, and we discuss the behaviour of the Lyapunov vector. 
Section [HI is devoted to a brief discussion of the travelling wave solutions. Finally, in 
Section[7|we present our conclusions and we discuss several proposals for future research. 
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2. Coupled logistic maps 

2.1. The logistic map 

The logistic map is a non-linear evolution equation acting on a continuous variable x 
taking values in the unit interval [0,1]. The evolution is defined by iterations over 
discrete time, n = 0, 1, . . .: 

i n = /(i n _i) = n n _i(l - i n _i) , (1) 

with the parameter r taking values < r < 4. The time series has very different 
behaviour depending on the value of r. For < r < 1 the iteration approaches the fixed 
point x* = 0. For 1 < r < 3 the asymptotic solution takes the finite value x*(r) = 1 — 1/ r 
for almost any initial condition. Beyond r = 3 the asymptotic solution bifurcates, x n 
oscillates between two values x\ and x%, and the solution has period 2. Increasing 
the value of r other bifurcations appear at sharp values, i.e. further period doubling 
takes place. Very complex dynamic behaviour arises from the relatively simple non- 
linear map (0) in the range r e [3.57, 1]: the map has bands of chaotic behaviour, i.e. 
different initial conditions exponentially diverge, intertwined with windows of periodic 
behaviour. Surprisingly enough, some exact solutions are known for special values of 
r jHIj. For example x n = sin 2 (6^2™) with k integer and 6 determined by the initial 
condition through Xq = sin 2 (07r) is a solution for r = 4 [3*2"] . 

2.2. Coupled map lattices 

Coupled map lattices (CMLs) jH] are discrete arrays of scalar variables taking continuous 
values, typically in the unit interval, that evolve over discrete time according to a 
dynamic rule. They are generalisations of cellular automata for which the variables take 
only discrete values. Kaneko introduced them as phenomenological models to describe 
media with high energy pumping but they may also arise as discrete versions of partial 
differential equations. 

A typical realization of the evolution of a single independent unit is given by the 
logistic map (0). One usually uses models in one or two space dimensions with periodic 
boundary conditions (a ring or a torus). A Laplacian coupling among nearest-neighbours 
on the lattice is usually chosen as it is motivated by the intent to model fluid mechanics 
in a simpler manner. In d = 1, and labelling with i — 1, . . . , N the sites on the ring, the 
coupled map model takes the form 

< = /(4-i) + \ - 2/(4-i) + (2) 

with x l ^~ N = x l n for all n with the number of elements on the ring. The initial 
condition is usually chosen to be random and thus taken from the uniform distribution 
on the interval [0, 1] independently on each site, v is the coupling strength between the 
nodes and plays the role of a viscosity. 

Numerical simulations showed that CMLs exhibit a large variety of space-time 
patterns: kink-antikink configurations, space-time periodic structures, wavelike patterns 
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or spatially periodic structures with steady, periodic or chaotic dynamics, space-time 
intermittent and spatio-temporal chaos have been found for different values of the 
parameters (r and v). A detailed description of the phase diagram is given in [Sj. In 
a nutshell, the dynamics is characterised by a competition between the diffusion term, 
that tends to produce an homogeneous behaviour in space, and the chaotic motion of 
each unit, that favours spatial inhomogeneous behaviour due to the high sensitivity to 
the initial conditions. 



3. Partial differential equations 

The d — 1 Kuramoto-Sivashinsky (KS) equation was introduced by Kuramoto to study 
local phase turbulence in cyclic chemical reactions [22] . Ad = 2 extension of it was later 
used by Sivashinsky to study the propagation of flame fronts in mild combustion [23J. 
In d — 1 it reads 

dh _ d 2 h d 4 h ^dh 

dt dx 2 dx 4 dx 

with h = h(x, t) a real function and k a real parameter. Similar pattern formation to the 
one developed in CMLs has been observed in the numerical solutions of the KS partial 
differential equation 33 . 

The Kardar-Parisi- Zhang (KPZ) equation was proposed as a non-linear model that 
describes surface growth [21]. If h(x,t) is the height of a surface on a substrate point x 
at time t, the equation reads 

dh 8 2 h A (dh\ 2 , , 

with r](x, t) a Gaussian white noise with zero mean and correlation given by 

( rj(x, t)r](x', t')) = D 8{x - x) 5(t - t') . (5) 

The variable dh/dx satisfies a noisy Burgers equation |25j . (For a discussion of the 
properties of these equations see ref. 

Using a dynamic renormalisation group calculation, Yakhot suggested that the 
elimination of large wave-vector modes generates a random 'stirring force' with zero 
mean and average {fi(k,u>)fj(k',u')) = kikj 5(k + k') 5{u + lu'). The average (. . .) 
represents a time-average in the stationary state. Simplifying further the propagator in 
the k — > limit, he argued that in the large scale and long time limit, k — > and u — > 0, 
the KS equation becomes the random-force-driven Burgers equation or, under a change 
of variables, the KPZ equation (with positive viscosity in d = 1 and negative viscosity 
in d> 2) dU. 

Next came a number of papers in which L'vov, Procaccia et al ^H]-^! studied this 
problem in more detail. They showed that for both equations the field-field correlation 
in Fourier space, n(k,uj), defined from 

( h(k, u)h*{k', u') ) = n{k, u) 5{k - k') 5{u - u') , (6) 
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has a scale invariant form 

with the 'roughness exponent' a and the 'dynamic exponent' z being equal to 1/2 and 
3/2 in d = 1, respectively [T^j. Under the assumption that a scale invariant solution 
of the form ((7j) also exists in d > 1, L'vov and Procaccia ^Hj argued that the scaling 
solutions for the two models bifurcate in d — 2; that is to say, the exponents z and a 
are not the same in all dimensions and hence the two equations do not belong to the 
same universality class. Further evidence for the breakdown of the relation between KS 
and the strong coupling phase of KPZ in d > 2 appeared in [T7j . 

It has been difficult to verify the coincidence of the scaling solution to the KS and 
KPZ equations in d — 1 by solving the KS equation numerically and comparing to the 
prediction z = 3/2. The reason for this difficulty is the existence of a very long crossover 
regime [T%] . A careful numerical study appeared in jTHj where large scale simulations 
were confronted to the predictions for the crossover behaviour obtained from the analysis 
of the KPZ equation. 

The Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) non-linear diffusion equation 
determines the evolution of the concentration of some chemical species or individuals, 
< h < 1, on a one dimensional space and reads 
dh d 2 h „.,. 

with f{h) a non-linear term, typically of the form f(h) = k\h — k 2 h 2 with k\ > and 
k 2 > 0. The FKPP equation has travelling wave solutions [3U| 

h(x,t) = F(x-vt) (9) 

that represent the invasion of the stable phase h(—oo, 0) = ki/k 2 in the unstable phase 
h(oo, 0) = and travel with velocity v > 2yfk[. Which velocity is selected depends 
on the initial condition. In many cases, if the initial front is "sufficiently steep", i.e. 
F(y) decays faster than e~ lminV for y — > oo f the front advances asymptotically with the 
minimal velocity u m $ n . These are called "pulled fronts": the leading edge of the front 
pulls the interface through growing linear perturbations about the unstable h = value. 
In these cases, the velocity approaches v m i n as a power law 

v{t) = v mm - -1— + 0(t- 3 / 2 ) . (10) 

In "pushed fronts" instead it is the non-linear growth in the region behind the leading 
edge that pushes the interface and the asymptotic velocity is larger than v m i n - 



4. The continuum limit of the CMLs 



In this Section we define the continuum limit and we briefly discuss each term in the 
resulting partial differential equation by making an explicit comparison to the ones in 
the KPZ equation (@J) and the FKPP equation ©• 
| a step function realizes this requirement 
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4-1. The continuum limit 

The main idea is to take the continuum limit of the CML using the following 
discretisation of time and space derivatives: 



3h h l n+1 - hi 



dt St 
dh hl +1 - hi 



(11) 
(12) 



dx 8x 
d^h K+ 1 - 2h\ + h^ 

dx* (s x y ' 

&h - + m n - 4/c 1 + K 2 

dx A ^ (5x) A ' 1 ' 

with St the time-step and Sx the lattice spacing equal one in our system of units. The 
CML of logistic elements then becomes 

where we called x the coordinate (iSx — > x), t the time (nSt — > t), and /i the field 
[x l n — > h(x,t) = h]. 

4-2. Relation with KPZ 

Comparing to eq. (jlj) one notices that: 

(i) By definition the field h is bounded and takes values in the unit interval. Thus, the 
resulting equation should have an effective confining potential that limits the field 
to a finite range. 

(ii) The elastic term is here multiplied by a field-dependent viscosity 

u{h) = ^(1 - 2h) . (16) 

(iii) The second, non-linear term is of the form of the one in the KPZ equation with a 
negative coupling 

A = -vr . (17) 

(iv) The last two terms read 

r}(x, t) = (r- l)h(x, t) - rh 2 (x, t) . (18) 

We notice that these terms are not present in the KPZ equation. In order to 
compare to the latter we shall argue that they have a double identity: on the one 
hand r] behaves roughly as a short-range correlated noise in space and time; on the 
other hand it can be interpreted as a force derived from a confining potential 

, = -^, ™ = -^» + ^. (19) 

In the next subsections we briefly discuss the elastic and non-linear terms. The 
analysis of 77 is more delicate and we postpone it to the next section where we present 
the results of the numerical integration of the CML. 
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4-2.1. The elastic term. The fact that the elastic term is multiplied by a field- 
dependent viscosity needs a careful inspection. 

First, v(h) is negative for h < 1/2 which implies an instability in the hydrodynamic 
limit. The same feature appears in the KS equation [eq. (jSJ)]. It was shown that this 
instability taps the system and so creates an effective 'noise' leading to the mapping of 
the KS equation onto a similar equation with a noise term and a renormalised positive 
viscosity coefficient [TH [T3] (see the discussion in Sect.EJ). With respect to this 'noise 
creation' feature, it is also very important to have a confining potential which restraints 
the instabilities caused by the negative values of v(h). Together they create the effective 
noise. Therefore, the fact that we may get negative values for the viscosity actually 
supports our line of argumentation, i.e. that of relating a deterministic model to a 
stochastic one. 

Second, if h remains bounded between, say, and 1 the 'bare' viscosity takes values 
on the finite interval [-y, y]. However, as is well known, in the KPZ system the 
viscosity coefficient is renormalised by the nonlinear term (put in simple words, the 
nonlinear term has a smoothing effect) so that the large-scale viscosity that the system 
experiences is not only determined by the bare value. Thus, we also expect the field- 
dependent bare viscosity to be renormalised at large scales and thus its precise value 
not to be very important. 

4-2.2. The non-linear term. It is well-known that the sign of the coupling constant is 
not important in the KPZ equation. The reason is that an inversion of the sign of A 
just corresponds to describing the evolution of the mirror image of the original surface 
(i.e. a description of —h(x,t), see e.g. [26J. Actually, when mapping the noisy Burgers 
equation onto the KPZ equation the resulting coupling constant is negative, namely 
A = -1. 

4.3. Relation with FKPP 
Comparing to eq. (jSJ) one realizes that 

(i) The field h is bounded as in FKPP. 

(ii) The viscosity is now field dependent and may take negative values. 

(iii) There is a KPZ-like non-linearity, not present in the FKPP equation. 

(iv) The last two terms are identical to the ones in the FKPP equation with k\ = r — 1 
and k 2 = r. 

The similarity with the FKPP equation suggests to search for travelling wave 
solutions in the CML and study their properties (velocity, etc.) as done for the FKPP 
equation (HHj. We shall come back to this issue in Sect. H3 
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5. Numerical tests in the context of surface growth 

We integrated numerically the CML of logistic units with up to N = 1024 sites and 
periodic boundary conditions. All units were updated in parallel. We used floating 
point precision 10~ 16 . Unless otherwise stated we chose independent random initial 
conditions for each map taken from the unit interval with flat probability. Since finite- 
size effects seem to be negligible for systems larger than approximately 32 sites [45] . we 
do not discuss smaller systems, nor do we make a systematic study of finite-size scaling. 
In particular, we discuss the choice of parameters in Sect. 15.11 we study the statistical 
properties of the r] term in Sect. 15 .2\ we compare the spatial and temporal behaviour 
of correlation functions in the CML and KPZ growth in Sect. 15.31 and we analyse the 
Lyapunov vector in Sect. 15.41 

5.1. Choice of parameters 

The panels in Fig. ^present the space-time plot of a coarse-grained variable obtained by 
transforming the continuous variable x l n into a bi- valued Ising-like one s l n = sign(x^— x*) 
for several values of the non-linear parameter r and the coupling strength v = 0.4. We 
draw time on the horizontal axis and space on the vertical one. Every 50 time step is 
plotted; if x\ is larger than x* = 1 — 1/r (the unstable fixed point of a single logistic map) 
the corresponding pixel is painted black; otherwise it is left white. This analysis allows 
us to identify different regions of phase space [H] in which we can study the connection 
with the KPZ equation. In this work we focus on v = 0.4 and r = 4, i.e. deep in the 
chaotic regime. We briefly mention at the end of this Section the behaviour found for 
other values of v and r. 




2500 5000 7500 10000 2500 5000 7500 10000 2500 5000 7500 10000 

Time Time Time 



Figure 1. Space-time plots of the coarse-grained variable s % n — sign(x^ — x*) for (a) 
r = 4, (b) r = 3.82, and (c) r = 3.8. 

5.2. Analysis of r)(x,t) 

5.2.1. The noisy aspect. Before studying the properties of t] in the CML let us recall 
some properties of this quantity in the single map. Choosing the parameter r to be r = 4, 
i.e. well in the chaotic regime, we obtain the histogram of x n over a time window of 
length T = 10 6 shown in Fig.Efa). The figure demonstrates that the numerical precision 
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chosen is good enough for our purposes since the data are rather well described by the 
analytic prediction for the probability distribution function (PDF) 

1 



p(x) 



7T\ x(l — x) 



(20) 



The divergence of the peaks at x = and x = 1 is however suppressed due to the floating 
point precision of the numerical data. The histogram is symmetric around x = 1/2 and, 
consequently, the mean is given by (x) = T _1 X)n-i x « = V^- 
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Figure 2. (a) The PDF of a single map x n , and (b) the evolution of its autocorrelation 
function defined in cq. (|2TJ> . 



In Fig. |2(b) we show the connected correlation over time, 

i T-n 

C x n = ( Ax k+n Ax k ) = V Ax k+n Ax k , (21) 

T - n kT 

for a single run using a randomly chosen seed. Ax/c = x k — 1/2 and the average is taken 
over all pairs of data obtained with a time delay n on a single run of length T = 10 6 . 
No correlations can be identified from this plot though we know, however, that some 
time-structure might be present [US] (see, e.g. the return map in Fig. |HJ). 
We now turn to the study of 

7] n = (r - l)x n - rxl . (22) 

By definition, r] n takes values in the interval [—1, (r — l) 2 / (4r)] that for r = 4 becomes 
[—1,9/16]. In Fig. 0(a) we show the histogram of r\ n corresponding to the same data 
shown in Fig. El Actually, using the PDF of x n in eq. (}2T)j) we derive the PDF of r] n that 
is given by 

1 
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where 9(rj) is the Heaviside (step) function. The numerical histogram is well-described 
by this analytic function with the proviso that the divergences are replaced by the peaks 
in the figure. Using this expression for p{rj) one finds that the average of 77 vanishes 




Figure 3. (a) The PDF of the 'noise' term in a single map r\ n . (b) The evolution of 
its connected autocorrelation function, C^, over time. 



1 3 

(77) = (r -l)(x) -r(x 2 } = 3x - - 4 x - = . (24) 

2 8 

A vanishing mean is recovered by computing the average numerically. Similarly we get 

(n 3 ) 3 (ri A ) 9 

Skewness = . = -- , Kurtosis = — — - = - . (25) 
(t7 2 ) 3 / 2 4 (r] 2 ) 2 4 

The time dependence of the connected autocorrelations, C% = ( Ar/k+n^Vk ) = 
(T - n)~ Efc=o Ar] k+n Ar] k =(T - n)' 1 J2k=oVk+nVk, are shown in Fig. EJb) and again 
no obvious structure is observed. 

Next, we study the noise term in the CML. In Fig. 01a) we show a representative 
histogram of the individual units x\. This PDF is clearly smoother than the one of the 
independent map shown in Fig. |2fa). For the rather high value of the coupling strength 
used, v = 0.4, all units behave statistically in the same way. The average is ( x ) ~ 0.67. 
In Fig. Efb) we show a representative histogram of the individual noise terms. Apart 
from a high peak at 77 ~ 0.5625 (which is just the maximum allowed value 9/16 within 
our numerical accuracy) the plot is almost flat. The average is (77) ~ 0. 

Since the connected correlation over time of a single map is rather noisy we prefer 
to show its average over all sites in the system. More precisely, in Fig. El we show the 
time decay of 

l N / x T-n \ 

[CJ] = [<AW^ fc )] = -5: — — £« + J , (26) 



where we used (rjk) ~ 0. Here and in what follows we use square brackets to denote an 
average over space. We chose the scaling with n 2//3 of the time axis and a logarithmic 
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Figure 4. The PDF of (a) a typical map x\ and (b) a typical noise term in the fully 
chaotic regime (r = 4). 



scale of the vertical axis to make clear that the time- correlations decay as a stretched 
exponential: 

[Cn] ~e" an2/3 . (27) 




5 10 15 20 25 5 10 15 20 



Figure 5. Decay in time of the auto-correlation of the noise on a site averaged over all 
sites in the sample, [C%], as a function of n 2//3 [(a) linear plot, (b) linear-logarithmic 
scale] . 

Even though we cannot exclude an additional power law dependence on time this result 
demonstrates that the time correlations of the noise r\ are very short-ranged. Note the 
oscillations at short times in Fig.|5^a) that are not visible in Fig.E^b) due to the absolute 
value. 

We also studied the spatial correlations of rj. In Fig. |Bl we display the equal time 
correlation 
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Here and in what follows the overline indicates an average over time. The line is a 
guide-to-the-eye showing the exponential decay 



f ~ 46 , (29) 

defining a 'correlation length', £, which is finite but rather long. The oscillations can 
also be taken into account, but we omit their precise functional description here. 




Figure 6. Decay of the equal-time spatial correlations - C] 7 denned in eq. (|28|l (in 
linear-logarithmic scale). The line corresponds to e~ 1 ^ with £ w 46. 



The short-range correlations in time and space shown in Figs. Eland El are confirmed 
by an analysis of the PDFs of space and time averages of the noise itself, 



1 



7? = ^ £ Vn 



. Vr, 



n=l 



1 N . 

T7 X] ' 

iv i=l 



(30) 



that are quite Gaussian, as shown in Fig. Ufa) and Fig. |7fb), respectively. The data 
used to draw the histograms are taken from the horizontal average at fixed vertical 
position, and the vertical average at fixed horizontal position, respectively, of the original 
continuous value of x\ that gave rise to the space-time map of Fig. d^a). One can 
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Figure 7. The PDFs of (a) 7f, (b) [r/„], and (c) Re [»$] for jfc = 10, 20 
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also draw the histograms of different Fourier components rf [the case k = actually 
corresponds to Fig. Efa)] which turn out to be also Gaussian, as can be seen in FigEfc) 
where the cases k = 10 and = 20 are displayed. 
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5.2.2. The confining aspect. In the previous subsection we studied the statistical 
properties of the values taken by rf n along a dynamical run. We here analyse the effect 
of deterministic force deriving from the potential (|19jh 

The fact that rf n is not really a perfect noise can be seen from the study of the 
return maps. Figure Efa) shows the first return map, i.e. the plot of rf n+x as a function 
of rf n for many values of n and a chosen site on the lattice. The information encoded 
in this plot is not included in the analysis of the temporal correlations. Even though 




Figure 8. (a) First return map, rf n+l against rf n , for a typical site in the fully chaotic 
regime (r = 4). (b) The confining potential V(h) given in eq. (|19J) . 

the map is not a simple one-dimensional line, but rather occupies some non-vanishing 
area, it does not fill phase space. In other words, the return map is not structureless 
and therefore in some sense not completely chaotic/stochastic (see ref. This plot 

demonstrates that unpredictability (i.e. no correlations) does not imply randomness. 

The bound on the original logistic map variables should translate onto a confinement 
of the surface fluctuations h. The interpretation of 1] as a force deriving from the 
potential V(h) shown in Fig. IHfb) allows one to identify this term as an important one 
providing the confinement of the fluctuations. For small h, more precisely, h < (r — l)/r 
(=3/4 for r = 4) the force is positive and it tends to increase the value of h. Instead, 
for large h, h > (r — l)/r the force is negative and it pushes the fluctuation to take 
smaller values [a similar confining effect is produced in the KS equation @ by the term 
—kK\. 

5.2.3. The non-linear aspect. 77 has the same structure as the non-linear terms in the 
diffusive and otherwise linear FKPP equation. It is then clear that eq. ()15|) admits the 
same spatially uniform and constant fixed points h = and h = (r — l)/r. One can then 
expect to find travelling wave solutions when special initial conditions are chosen. We 
present some examples in Sect. EH but we delay their detailed numerical study to [56] . 
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5.3. The CML against KPZ: observables 

We here translate the definition of the observables of interest in the KPZ framework to 
that of the CML and we confront their behaviour in the two system. 

5.3.1. Confinement. In the CML the variables are bounded and the mean-square 
displacement cannot be larger than 1. In Fig.|2]we show a snapshot of a CML (a), a KPZ 
surface (b) and a confined KPZ surface (c). At face value these figures look different, 
although one can find similarities between the CML and confined KPZ surfaces. Also, 
these figures are very different from the KPZ surface being confined between two 'soft' 
walls 36J . 




(a) J 

20 



„ 




Figure 9. (a) A CML interface, (b) a KPZ interface and (c) a confined KPZ interface 
(with a confining potential of —-^h 2 ). 

Let us now compare the space and time dependence of these surfaces in more detail. 



5.3.2. Roughness. Stochastic surface growth from a flat initial configuration is usually 
characterised in terms of the mean-square displacement w 2 . In a one dimensional 
discrete model w 2 reads 

2 



h\t) 



L 3=1 



(31) 



( w 2 {L, t) ) has the 



u <C 1 (growth) , 
(saturation) 



(32) 



The square root of the noise averaged displacement, w(L,t) 
scaling behaviour j2E] 

w{L,t) ~ L a f (—) , f(u)~\ UZ ' 
y ' J \L Z ) ' 1 const , u > 1 

with a and z the roughness and dynamic exponents, respectively. 

Since the CML surface is necessarily bounded, if the scaling form holds, a 
must be identically zero. If we wish to check z numerically we need to specify the initial 
conditions. In the CML a flat configuration for all sites remains the same at 

all later times. In order to generate a rough configuration we need to introduce some 
disorder initially. This can be done in a number of ways. For example, one can analyse 
the motion of an initial localised bump on an otherwise flat configuration or one can add 
some small random noise on a flat state. In Fig. El we present results of a simulation 
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of the CML, with random, though very small, initial condition and for several system 
sizes. 



-1 r 




1 2 3,4 5 6 7 

In n 



Figure 10. Roughening of a CML with small random initial conditions for various 
system sizes given in the key. 

As can be seen in the figure, the saturation time does not depend on the size of the 
system, and occurs after ~ 10 time steps in all cases. Also, the width in the steady-state 
(i.e. w n for large n's) is also independent of system size. These findings are consistent 
with vanishing exponents z = and a = [see eq. ()32)) ]. Actually, vanishing a and z is 
also consistent with a simple scaling argument applied to the confined KPZ equation, 
with a confining potential of — |/i 2 (at this point it is worthwhile mentioning that some 
confined KPZ systems have been studied in the past jHHj though, as far as we know, 
non is directly relevant here). 

The non-trivial character of the PDF of mean-square displacements found in several 
unbounded surface growth problems [37] is completely erased by the confinement. 



5.3.3. Stretched exponential relaxation. The decay in time of the field-field correlations 
is still non-trivial in this system and very similar to the one found for the (usual and 
unbounded) KPZ equation. In Fig. [TT] we show the space averaged, time correlation 
of the field [ ( h m+n h m )]. The linear dependence in the linear-logarithmic scale used in 
panel (b) indicates the stretched exponential decay 

tt(t)^e-" C , with C = 2/3, (33) 

The panel (a) shows that the decay occurs in an oscillatory way. 

A stretched exponential relaxation of the dynamical structure factor in the d — 1- 
KPZ growth, 

C{k, t) = ( h{k, t)h*{k, 0) ) ~ Ae- Bktl/Z (34) 

with z the dynamic exponent was found by solving the KPZ equation within the mode- 
coupling approximation [HE], with a self-consistent expansion |39J and with a direct 
numerical integration in d = 1 jJJ]]. An oscillatory superimposed dependence on time 
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was obtained numerically |3*H1 |4TJ] and can also be obtained analytically from the self- 
consistent expansion [H]. Some doubt on this result was shed by the solution to a 
polynuclear growth problem that is in the same universality class as d — 1 KPZ [12] 
with a dynamic structure factor with a simple exponential decay. This discrepancy 
may be due to the fact that defining universality classes out of equilibrium might 
tricky: some models may share the static exponents but differ in their overall dynamical 
behaviour |%H] . 

We thus found that the dynamic structure factor of the CML of logistic units 
in their deep chaotic regime (r = 4) and the d = 1 KPZ equation have the same 
stretched exponential relaxation. Surprisingly enough, the dynamic exponent z = 3/2 
is reproduced in the CML model where it does not have the interpretation of relating 
spatial and time fluctuations: the study of the mean-square fluctuations yields a trivial 
size independent result (z = 0) due to confinement. The dependence of the exponent £ 
on the non-linear parameter r is weak, as demonstrated in Fig. IT2l where we display the 
decay of the dynamic structure factor at k = for several values of r. 




10 20 0499 30 40 50 10 2 °0611 30 40 50 5 10 ^jj^ 20 25 30 



Figure 12. The r dependence of the stretched exponential relation for (a) r = 3.83, 
(b) r = 3.87 and (c) r = 3.9. 

In a sense, it is not surprising that the bound on the values of the fluctuating field 
does not affect the time decay of the correlations. In the context of glassy systems 
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stretched exponential relaxations of correlation functions were searched and found in 
kinetically constrained models m which the dynamic variables are naturally bounded 
(they are spins). 

It would be very interesting to test whether a similar stretched exponential 
relaxation occurs in the confined continuous KPZ equation but this check goes beyond 
the scope of this article. 



5-4- The Lyapunov vector 

It was conjectured by Pikovsky et al [21] that the evolution of linear perturbations 
in the chaotic regime of a CML of logistic elements is well described by the KPZ 
stochastic dynamics. The notion of a Lyapunov vector is one of the ways to extend 
the notion of Lyapunov exponent to space-time chaos. One solves simultaneously the 
original nonlinear dynamical system, and the linearized equation for a perturbation. 
The Lyapunov exponent is then determined from the norm of the Lyapunov vector. 

In the case of the CML, it is convenient to write the linearisation of eq. (J2J in the 
following way 

u( x , t) = f [h(x, t)] u(x, i) , (35) 
u(x, t + 1) = (1 - v) uj{x, t) + - [Q{x -l,t) + uj{x + 1, t)} , (36) 

where u(x,t) is the Lyapunov vector, using the continuum limit (fTTj) - (fTlj) we obtain 
the following continuum equation for Cu(x,t): 



8Cj v d 2 u , „, r , . 

+ {/' [h(x,t + 1)]-1} 



(37) 



dt 2 dx 2 

We are then led to identify a 'multiplicative noise' 

Z x (x,t)=f'[h(x,t + 1)]-1. (38) 

(The fact that /' [h (x,t + 1)] is strictly speaking evaluated at time t + 1 is not a real 
problem since the process does not depend on the evolution of uj(x,t), as it does not 
feed into the equation for the h's. Therefore, one can say that intrinsically the noise in 
our system should be understood in the Ito sense.) 

At this point, since the variable u(x,t) is not bounded, a simple scaling argument 
can convince us that the coupling of the second derivative to the noise £ x is less relevant 
(in the RG sense) than the coupling to uj(x, t) itself. Therefore, we neglect the last term 
inside the square brackets. The equation we are left with is just a diffusion equation 
in the presence of multiplicative noise. This motivates the application of the Hopf-Cole 
transformation, H = ma), which yields 



dH ud 2 H v (dH\ 



2 



m -2 3* • 2{T*) +^ t >' < 39) 

Now, we have to figure out the dynamics of H = In a;, i.e. the Lyapunov vector. 
This is easily done by taking the In of eq. (J35j) 

H(x,t) = H(x,t) -\nf[h(x,t)\ (40) 
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This means that H will be just H up to an additional 'additive noise', (x, t), defined 

as 

Z + (x,t) = -]nf[h{x,t)} . (41) 

The roughness of the 'interface' H(x,t) is given by the equal-time autocorrelation 
function of H(x,t). Therefore, we need to check the statistical properties of £ + (x,t), 
that is to say, its autocorrelation function and its correlations with H(x,t). We checked 
that since the term is exponentially distributed, with very short-range correlations in 
space-time and vanishing correlations with H(x,t), the dominant contribution to the 
autocorrelation function of H(x,t) comes from H(x,t), and is therefore described by 
the KPZ exponents. This is not surprising, since the noise £ + (x,t) does not accumulate 
in H(x,t), but rather it just adds a random number to H(x,t) with a small amplitude. 

This establishes the hypothesis of Pikovsky et al |2Zj but, more importantly, it 
allows us to make a step forward relating known quantities in the KPZ literature to 
features of the Lyapunov exponent. The Lyapunov exponent is obtained from the norm 
of the Lyapunov vector. If one uses the so-called 0-norm, 

1 r L 



N (t) = exp 



h(x, t)dx 



(42) 



L Jo 

it is then assured to be a self-averaging quantity. In addition, the Lyapunov exponent 
is given by 

hAMD-hAMO) 

T^oo T 

and this is no other than the large-deviation function for the Asymmetric Exclusion 
Process (ASEP). calculated previously by Derrida and Appert [13]. The ASEP is a 
discrete model in the universality class as KPZ in one dimension. In terms of the ASEP, 
the Lyapunov exponent is given by 



'p(l-p). 



A = p (1 - p) + y ^jjrG [y]2*p (1 - p) Lj (44) 

where L is the system size, p is a the density of particles (a parameter in ASEP), and 
G(0) is a scaling function independent of L and p and known is an implicit form (see [13] 
for more details). 

What is left now is to relate the parameters of the ASEP model to those of the the 
KPZ equation we derived for H(x,t) (JSDJ) - Using ref. [Sj we see that p(l — p) = D/32 
where D is the noise amplitude characterizing £ x (x,t). In general, D would be a 
complicated function of r and u, and is unknown at the moment. Still, the expression for 
the Lyapunov exponent ([Hjl could be useful to study finite size effects. As an exercise, 
we try to estimate D using the known PDF of a single (uncoupled) map, in the fully 
chaotic regime (r = 4), given by eq. ([27!|) . In that limit £x = 3 — 8x, and is known to 
be uncorrelated in time. Therefore, D = = 9. Plugging this into eq. ([4^[) yields 
A = 9/32 ~ 0.28. In numerical simulations Pikovsky and Politi found A ~ 0.38 [2Z] 
which is of the same order of magnitude as our result. The difference is certainly related 
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to the fact that we approximated £ x by that of a single map. In addition, we neglected 
£+, and this could have influenced as well. 

6. Travelling wave solutions 

Let us now briefly analyse the coupled map lattice model and its continuum limit from 
the viewpoint of front propagation and look for travelling wave configurations linking 
the two stationary solutions: 

h = , and h= T -^- . (45) 

r 

The non-linear term rj renders the solution h = linearly unstable: dr]/dh(h = 0) = 
(r — 1) > and the solution h — (r — l)/r stable drj/dh(h = (r — l)/r) < if r > 1/2. 
Introducing the travelling wave Ansatz (jOJ) in the linearized version of eq. in which 
we drop the terms proportional to F 2 , FF" and {F') 2 we find 

v > v min with v min = \Jlvr{r - 1) , (46) 

and 7 m i n = ^2(r — l)/(z/r). The combined effect of the diffusion and non-linear term 
of KPZ type are not obvious but one may expect that if we start with a localized 
initial condition with, e.g., the form of a bump with h ^ in a finite region of the 
one dimensional space, the non-linear term will drive the border of the bump into the 
unstable state h = 0. The front may then advance as a travelling wave. 




Figure 13. Evolution of the initial configurations for (a) a step-like initial condition 
and (b) a Gaussian bump. 

In Fig. EI] (a) and (b) we show the evolution of the CML starting from localized 
configurations with a step-like initial condition in panel (a) and a Gaussian form in panel 
(b). The configuration x l n at three times, n = 0, n = 200 and n = 400 are shown with 
different lines. They have been translated vertically to render the identification of the 
curves simpler. For both initial conditions, after a very short transient the configuration 
acquires rather sharp borders with a 'random'-like form within the step. More precisely, 
the value at the plateau close to the borders is well described by the stable solution 
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x n ~ (r — l)/r = 3/4 while the values of the random-looking peaks on this plateau have 
an average ( x ) ~ 0.67 as in Sect. 15.21 For these times the decay to zero at the edges is 
very sharp, it occurs in less than 10 lattice spacings. 

In Fig. E]we display the evolution of the CML starting now from a more extended 
initial conditions; more precisely, we took an exponentially decaying form e~ x ^ with 
£ = 8, and we show the configuration x l n at three times n = 0, n = 20, and n = 40. It 
is clear that the front travels much faster in this case than in the examples studied in 
Fig. E2 Moreover, the interface is not as sharp. 




Figure 14. Evolution of the initial configurations for an exponentially decaying initial 
condition. 

We delay for future work [Slj a quantitative analysis of the velocity of propagation. 
One might expect that additional non-linearities, as the KPZ-like term, can only increase 
(or leave unchanged) the velocity with respect to v min . A sufficient condition to have 
a pulled front is that all non-linear terms should suppress growth for a front that 
propagates into a linearly unstable state. In our modified FKPP equation it is not 
obvious a priori whether this holds. 

7. Conclusions 

We have shown that in the continuum space and time limits the CML of logistic units 
in d = 1 becomes an equation of the KPZ type. We argued that some non-linear 
terms have a triple effect. These terms can be interpreted as a noise source with 
short-range correlations in time and space, an assumption that we supported with the 
statistical analysis of numerically generated configurations. But they also provide a 
source for confinement thus reducing considerably the roughness of the surface: the 
positive roughness exponent a found in conventional KPZ is reduced to a vanishing 
value. On the other hand, the same non-linear terms have the form encountered in 
the FKPP equation for the diffusion of advantageous genes and one can then expect 
to find travelling wave configurations in the discrete model. We also argued that the 
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field- dependent viscosity appearing in the partial differential equation has no special 
effect. 

As regards to time-dependent observables in the CML, we found that the 
correlations in time of the local fluctuations have a stretched exponential decay as 
recently claimed to arise for surfaces generated by the conventional (unbounded) KPZ 
equation |3E|-|I3- 

We have shown results for the pair of parameters v — 0.4 and r = 4 where the 
system is deep in the chaotic regime. A detailed analysis of the phase diagram is 
beyond the scope of this article. We just mention that the stretched exponential decay 
persists when keeping the value of v fixed and reducing r until reaching the 'critical' 
value r « 3.83. Interestingly enough, long periods of local blocking are not necessary 
to find such a slow relaxation since we find it even when the space-time plots look 
extremely chaotic. Time correlations are hidden in these plots. 

A number of authors have signalled the possible relevance of CMLs in describing 
different aspects of glassy relaxation. Recently, Mousseau et al jlSl EE] studied the 
same model in its intermittent regime [y = 0.4, r w 3.83) with the aim of relating the 
10 orders of magnitude stretched exponential decay of its distribution of trapping times 
(times in which an element remains locked into one of the coarse-grained values s l n = ±1) 
to the one observed in super-cooled liquids. Simdyankin and Mousseau [IE] associated 
this stretched exponential decay to the one of the correlation function. As we have 
already stressed, we also obtain a stretched exponential relaxation for larger value of 
r where trapping intervals for the coarse-grained variables s l n do not exist (see Fig. QJ. 
In a similar spirit to the works of Mousseau, Garrahan and Chandler associated jSU] 
the slow dynamics and glass transition to the structure of trajectories, in the form of 
the space-time maps, of the spin-like variables in kinetically constrained lattice gases 
(see jSI] for a review of these models). A fully-connected model of logistic maps was 
studied from the glassy point of view by a number of authors; a range of values of r with 
a large number of 'macroscopic states' and 'ergodicity breaking' jlTj and the possibility 
of an analog of replica-symmetry-breaking [IE] were reported. With the baggage gained 
from the current understanding of the dynamics of glassy systems jS2] we intend to 
revisit the generation of an effective temperature [H EH IHH] and its possible appearance 
in a fluctuation relation [EH EE] in chaotic systems [EE] • 

Finally, we also briefly analysed the resulting partial differential equation in 
comparison with the FKPP and we showed that the CML has travelling wave solutions 
with similar qualitative properties to the ones in FKPP. A more detailed analysis is 
necessary to understand the main characteristics of the front and, in particular, classify 
its velocity depending on the initial conditions and other parameters. 
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